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(54) Abstract Title 

Hybrid-linear-bicubic interpolation method and apparatus 

(57) The Invention relates to a method for reducing the computational overhead of bicubic interpolation while 
still providing a similar level of accuracy. The invention takes into account the fact that sampled points 
surrounding a point whose value is to be determined have respective first, second and third order effects on 
the calculated value. The Invention combines linear interpolation, Ignoring points having a third order effect, 
with cubic interpolation of points having a first and second order effect to derive the value. 
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1 

HYBRID -LINEAR -BICUBIC INTERPOLATION METHOD AND APPARATUS 

The present invention relates to a method and component for 
interpolating data. 

in the area of image and video processing, there is frequently a 
need to create new images by re -sampling an original image. 

Typical instances are: 

a) zooming- in, where there is a need to create more pixels than in 
the original image; 

b) image-rotation, where new pixels overlap 2 or more of the 
original pixels; and 

c) Image -distortion (morphing) , where the~aim may be to correct 
lens aberration (barrel, pin- cushion distortion) or even to 
deliberately introduce distortion. 

Several two-dimensional techniques for computing new pixel -values 
are well known, and used frequently. The most commonly met ones are: 

a) Nearest -neighbour 

b) Bilinear 

c) Bicubic 

Nearest -neighbour is extremely simple. It consists of sampling the 
pixel nearest to the centre of the desired new pixel, and using that 
value. Unfortunately, it has the undesirable effect of introducing 
'jaggies' or ' staircasing' . Consider the case of a horizontal line, 1 
pixel wide. This appears to the eye as a perfect, straight line (which it 
is). If re-sampled using nearest -neighbour , a very-slight image -rotation 
of a degree or so will convert this line into a highly visible staircase. 
Further operations will enhance this effect, and can even cause break-up 
of the original line. This form of interpolation is not considered 
further, as the quality is very poor. 

Bilinear is more complex, it consists of taking the 4 pixels which 
'surround' the position of the new pixel, and computing a new pixel -value 
based on these 4 neighbours, whichever pixel is nearest to the new one 
has most weight. If the new pixel is right in the middle of the 4 
neighbours, then it will end up being the average of the 4 pixels. This 
technique is a considerable improvement on nearest neighbour, but still 
has some draw-backs. It has the effect of blurring the original image, 
and losing contrast, as the new pixel's computed value will never lie 
outside the min-max range of the 4 neighbours. 
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Bicubic is recognised as being the 'best' quality, as it applies a 
curve-fit to the pixel-values, and so is attempting to re-construct the 
original waveform which gave rise to the sampled pixel -values in the 
original image. Unfortunately, bicubic interpolation is recognised as 
5 being expensive computationally. 

The present invention provides a new type of interpolation, a 
hybrid linear-bicubic interpolation, which achieves the quality of a 
bicubic interpolation, but using some 80% of the computation of 
10 traditional bicubic interpolation. 

According to the invention there is provided a method of 
determining an interpolated value for a point according to claim 1. 

An embodiment of the invention will now be described with reference 
to the accompanying drawings, in which: 

Figure 1 shows a waveform, which has been sampled at points WO, wl, 
W2, W3 to generate pixels in an original digital image; 

Figure 2 illustrates a point value v calculated according to linear 
interpolation; 

Figure 3 shows a 2-D grid of points; 

Figures 4 and 5 illustrate the location of a point V, whose value 
is to be calculated, within respective 2-D grids; 

Figure 6 illustrates the sequence of interpolations according to 
the invention to calculate the value at point V; and 

Figures 7 to 9 show the performance of the invention versus prior 
art methods in terms of Instructions Count, Memory- Fetch (load) 
instructions and LOT space required. 

in order to assess the present invention, we need to look at the 
three prior art techniques using straight calculation as well as 1- 
dimensional and 2-dimensional lookup-table (LUT) implementations which 
can be used in practical situations to speed-up the brute-force required 
by straight calculation. 

K little introduction is necessary here to show how one -dimensional 
linear and cubic interpolation actually work, with and without the use of 
lookup- tables, and to show the use of simple pseudo-code instructions 
(such as are found in microprocessors and DSP (Digital Signal Processing) 
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chips) . 



Consider the problem of interpolating new points between existing 
points, and, in particular, consider a one -dimensional case for 
simplicity. Figure 1 shows a waveform, which has been sampled at points 
WO, wi, W2, W3 to generate pixels in an original digital image. Note that 
the sampling rate is greater than twice the natural frequency, so the 
Nyquist requirement is observed, and the frequency can be captured 
without aliassing it. Because the original image has pixels of equal 
size, the sampled-values WO, wi, W2, W3 lie at equal intervals along the 
T-axis. This simplifies the maths, but is not a fundamental limitation on 
any of the interpolation methods described here. 

The requirement now is to generate new pixel values at points along 
the T axis, between the existing wi and w2 points. 

1-D LINEAR INTERPOLATION IN GENERAL 

The technique takes a straight line between wl and W2, and picks 
off as many points along it as are required. It is convenient to 
introduce a parameter, here called j, which is the distance along the T 
axis from the start of the arc (wl point) to the new required point, 
Figure 2. 

Also, let N be the distance along T axis between Wl and W2 points. 
N will in general be a known, pre • determined constant value. 

v is the new point whose value is to be computed. (The fact that v 
is a long way from the original curve is simply an indication of the poor 
quality, corner-cutting behaviour inherent in linear interpolation) . The 
origin of the T-axis is relocated so that: 

t=0 at w=wi, t=N at w=w2, t=j at w=v 

it is easy to show that: 

V « (j*W2 + (N-j)*Wl) / N 

Consider now the implementation of this technique with regard to real- 
world requirements, eg the availability or not of floating-point hardware 
and the use or not of lookup- tables (LUTs) . CPU registers will be called 
names like RO, Rl, R j , RN etc. 

N features as a known constant, which is assumed to be held in a 
register called RN which is read-only, as it is used so often. 
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EXACT 1-D LINEAR INTERPOLATION WITHOUT LOT 

Rj (=j) is assumed to be known, as a result of sqme previous 
morph/lookup/calculat ion . 

rn is known, wl and w2 are not yet known, and need to be loaded 
from memory in which the image is stored. Appendix #1 shows the pseudo- 
code to implement the technique. This uses 2 load-f rom-memory, 2 
multiply, 1 divide, 3 copy/add/sub instructions, total 8 instructions. 

in cases where further calculations are to be done, it's possible 
to work using N*v as a value, leaving it until later to factor-out N. 
This saves 1 divide instruction at this stage. Also, note that if it can 
be chosen that H-l. this also saves the divide instruction entirely. 

If the technique is implemented using floating-point hardware, than 
N =l 0 is easy to specify, if integer hardware is used, than N must be 
considerably greater than 1, as an integer value of j which lies between 
0 and n is required. In this case it makes sense to choose N to be some 
convenient power of 2, eg 8. 16. 256, whatever gives sufficient 
resolution. Using powers of 2 means the division by N can be achieved by 
bit-shifts (if cheaper than divides), or even avoided altogether by 
selecting high-byte sections of registers and ignoring low byte sections. 
This is architecture -dependent. 

in conclusion, exact 1-D linear interpolation without LUT requires 
8 instructions to calculate V, or 7 if N*v is wanted. 

1-D LINEAR INTERPOLATION USING 1-D LUT 

in the case of image pixels, typically integer values of w are 
used values from 0 to 255 are most commonly used, as in 8-bit grey-scale 
images and 24 -bit RGB and YUV colour images. The resolution of such 
values is one part in 256, so the use of floating-point arithmetic xs 
overkill in most cases, it is quite accurate enough to resolve the 3 
value to this same 1 in 256 resolution, in other words. N can be 256 and 
j restricted to have values 0, 1,2, . . . ,255. 

suppose j has been computed as j=127.567: the difference in 
interpolated pixel-values using floating-point maths and integer-maths 
with j truncated to 127 will not be visible. 

Thus, j is restricted to have a small, finite number of values 
(e o 0 to 255). j-fl will produce v-wl. j=255 will produce something very 
close to v=w2. (j=256 will produce v=w2. but this point should really be 
generated by interpolation between w2 and W3. with j-0) . 



For each value of j=0 to j=255, there is now a weighting that can 
be applied to Wl and W2 values, and a lookup- table can be used to get 
these, using j as the index into the LUT. As before, 

Rj (=j) is assumed to be known, 
RN is known (=N) , 

Wl and W2 are not yet known, and need to be loaded from the image memory. 

Appendix #2 shows the pseudo-code to implement the technique where 
Wl and W2 are multiplied by their weights, summed, and divided to get v. 
This requires 3 load-from-memory, 2 multiply, 1 divide, 1 add/sub 
instructions, in total 7 instructions. 

In cases where further calculations are to be done, it's possible 
to work using N*v as a value, leaving it until later to factor-out N. 
This saves 1 divide instruction at this stage. (Note: N cannot be chosen 
to be 1, as now using integer values, with N » 1) . 

In conclusion 1-D linear interpolation using a 1-D LUT requires 7 
instructions to calculate v, or 6 if N*v is wanted. 

Using a 1-D lut saves two arithmetic instructions, at the cost of 
one extra load. The size of the lut is of the order of N entries (denoted 
as 0(N) size) . 

1-D CUBIC INTERPOLATION IN GENERAL 

The technique takes a curve between wl and W2, and picks off as 
many points along it as are required. There are an infinite number of 
curves which can be drawn through wl and w2, so the problem is to choose 
a sensible one which is well-behaved, and amenable to easy computation. 
There are other requirements placed on the choice of curve; clearly, in 
the case of the linear interpolation above, there is a slope - 
discontinuity at Wl and W2 where the previous interpolated- lines join-up. 
In order to achieve smooth interpolation across new arcs (wO-Wl arc 
changeover to W1-W2 arc etc) it is required to have slope -continuity at 
Wl and W2. So, we have 4 conditions, slope -continuity at wl, W2, and 
position-continuity at wl and W2. The simplest algebraic curve that meets 
this is a. cubic curve. 

The choice of slope at wl and W2 is somewhat arbitrary; any rule 
for choosing must be symmetrical with respect to the previous arc, i.e. 
the slope chosen at the end of the previous arc must be the same as that 
at the start of the next arc. The simplest choice here is to choose the 
slope at Wl to be parallel to a line drawn between WO and W2, and the 
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slope at W2 to be parallel to a line between wi and W3. So, choose slope 
at Wl = (W2-W0)/(2*N), and similarly choose slope at w2 - <W3 -Wl> / (2*N) . 

Let "t**3" denote t- cubed. 
Let w=wi at t=0 
Let W=W2 at t=N 

Let "R2 R3" denote that R3 is added to contents of R2 register, 
using 'C notation. Similarly, "Rl /= R4» means Rl is divided by R4, etc. 
Let "R2 = LUT [Rj] " denote that Rj is used to index into a table of values 
(lookup table) called LUT, again using 'C notation. 

It can be shown that the section of cubic- curve between wi and w2 
with slopes defined as above is given by the equation: 

w = (A * <t**3>) + (B * (t**2)) + (C * t) + D 
where 

A = ( W3 - 3*W2 + 3*W1 - WO 1 / ( 2 * (N**3> ) 
B = ( -W3 + 4*W2 - 5*W1 + 2*W0 ) / ( 2 * (N**2) ) 
C = ( W2 - WO ) / ( 2 * N ) 
D = Wl 

BXACT 1-D CUBIC INTERPOLATION WITHOUT LUT 

As before: 

Rj (=j ) is assumed to be known, 
RN is known (=N) , 

WO. Wl, W2 and W3 are not yet known, and need to be loaded from the image 
memory. 

with simple arranging, it can be shown that we have: 

2*(N**3)*A = (W3-W0) - 3(W2-W1) 
2*(N**2)*B = (W2-W1) - 2(N**3)A - Wl + WO 
2*N*C = (W2-W0) 

Appendix #3 shows the pseudo-code to implement the technique. This 
requires 4 load- from -memory, 6 multiply, 1 divide, 13 copy/add/sub 
instructions, in total 24 instructions. 

Note that some multiplies could be replaced by shifts, if N is a 
power of 2, depending on CPU instruction- times and the slight loss in 
freedom to choose N to be any integer/floating-point value. 



If floating-point operations are used and n is chosen to be 1.0, 
then this eliminates 2 multiplies, reducing the count to 22 instructions. 

If 2(N**3)v can be used in later calculations, then a further 1 
division is saved. Using floating-point registers this is no problem. 
Using integer registers it's likely that this will be a large number, so 
requiring the division by 2<N**3) or at least by N**2 to reduce it to a 
manageable scale. 

In conclusion, exact 1-D cubic interpolation without LUT needs 24 
instructions to calculate v, or 23 if 2(N**3)v is wanted without needing 
any LUT -space. 

1-D CUBIC INTERPOLATION USING LUT 

As per the previous linear interpolation using a lut, n is 
specified to be a pre -determined fixed number (eg, but not 
necessarily=256) . For each value of j from 0 to N-l we can solve for 
A,B,C,D and substitute for j to get v. For example, at t = N/4 it may be 
shown that 

V (at t=N/4) = ( -9W0 + 111W1 + 29W2 -3W3 ) / 128 

In other words there are 4 weights which need to be multiplied by 

WO through to W4, summed, and divided to calculate v. (if floating-point 

is used, the weights can be pre-divided, saving the final division) . 4 

LUTs are needed, one each for WO, wi, W2, W3, and each is of size O(N). 
Call these LUTs LUTO, LUT1, LUT2, LUT 3 respectively. 

Appendix #4 shows the pseudo-code to implement the technique. This 
requires 8 load- from -memory, 4 multiply, 1 divide, 3 add instructions, in 
total 16 instructions. 

in conclusion, 1-D cubic interpolation using LUT (integer) needs 16 
instructions to calculate v, or 15 if 2(N**3)v is wanted and needs 0(4*N) 
entry LUT -space. 

2-D INTERPOLATION IN GENERAL 

Consider now the 2 -dimensional situation based on a 2-D grid of 
pixel -values. In particular, consider the bicubic situation, as the 
bilinear case is a simpler subset of the bicubic one. 4 rows of pixels 
are required, each comprising 4 values. Let these rows be called w, x, Y, 
z rows, with values W0-W4, X0-X4, Y0-Y4, Z0-Z4 respectively. 
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Figure 3 shows a 2-d grid of points. Imagine a flexible rubber- 
sheet surface passing through 4 points XI. X2, Yl. Y2. The value at each 
point on the surface corresponds to the height of the sheet at that 
point, we now need to interpolate one (or more) new points at some pre- 
calculated point on this surface, as in the 1-D case above, j is defined 
to be the distance along the Row-axis from the Wl/Xl/Yl/Zl column, and 
similarly k to be the distance down the column axis from the X0/X1/X2/X3 



row. 



Bicubic interpolation requires doing 4 cubic interpolations along 
each of the 4 rows using the same j value, to calculate 4 new 
intermediate points called P.Q.R.S in Figure 4. Using these 4 points one 
further cubic interpolation is done in the column- direction using the k 
value to get the final value V. 

Bilinear is similar, but simpler, as only points Q and R are 
calculated using linear interpolation between XI and X2, Yl and Y2, then 
q & R are used to calculate v. 

Turning to Bicubic interpolation, the use of LUTs presupposes the 

v1 vo Yl Y2 is divided into a grid, with a known N number 

and^ew pixels will be taken from nearest grid-points, as before. 

If LUTs are used to assist the calculations, then they can be used 
in two main ways: 

Firstly, in a 1-dimensional way, as for the earlier 1-D 

'■^mliations weights can be looked up for points NO, Wl, W2. W3. 
interpolations, weig game 

These same weights are used for rows X, Y, Z, as the vaiu 3 

L all 4 rows. New weights must then be looked up for points P.Q.R.S as 

l lil be different from So. using 1-D LUTs, there is a lot of re -use 

of the previously looked-up wei.hts. This requires 1 LUT of size of» 

entries. 

Secondly, in a 2-dimensional way, it can be shown <^°*>^ t0 
the 1-D case) that for each value j and k, there is a weight to be 
allied to each of the 16 points wO through to Z3. This 2-D lookup table 
neeas *"2 entries, and there are 16 such LUTs needed, one for each of XO 
To gh to Z3. (Symmetry can be used to reduce this to 3 LUTs > at *e 
cost of some index-manipulation at run-time) . This requires 16 LUTs of 
size 0(N**2) entries. 

The same applies to the simpler bilinear interpolation, but only 1 
LUT of size O(N) entries is required if using 3 1-D linear 



interpolations, or 4 luts of size 0{N**2) entries if using 1 single 2-D 
interpolation. 

Consider now the instructions required to perform 2-D interpolation 
using straightforward calculations, bilinear and bicubic methods, 1-D and 
2-D LUTS. 

EXACT BILINEAR INTERPOLATION WITHOUT LUT 

Referring to exact 1-D linear interpolation without LUT, and Figure 
5, 7 instructions are used to calculate value N*Q between points XI and 
X2. 

Appendix #5 shows the pseudo-code to implement the technique. After 
getting Q along the X1-X2 row, N*Q will not be huge, so it can be assumed 
no division by N is required until interpolating along Q to R. This 
requires 4 load- from -memory, 6 multiply, 1 divide, 7 copy/add/sub 
instructions, in total 18 instructions. 

An alternative is to work it out brute- force algebraically, where 
it may be shown that: 

V = UN-j) (N-k)Xl + j(N-k)X2 + (N-j)kYl + jkY2) / (N**2) 

This can be implemented as shown in Appendix #6, where N**2 is 
either a known constant, or already loaded into a register. 

In all this uses 4 load-from-memory, 6 multiply, 1 divide, 7 
copy/add/sub instructions, in total 18 instructions. It is also more 
suited to built-in matrix instructions designed into the instruction-set 
of some chips. 

in conclusion, exact bilinear interpolation without LUT needs 18 
instructions to calculate v and needs no LUT space. 

BILINEAR INTERPOLATION USING 1-D LUT 

Referring to linear interpolation using 1-D LUT, 6 instructions are 
required to calculate value N*Q between points XI and X2 . 

Appendix #7 shows the pseudo-code to implement the technique, in 
all this uses 6 load-from-memory, 6 multiply, 1 divide, 3 copy/add/sub 
instructions, in total 16 instructions, so no improvement. 

The alternative is to work it out brute- force algebraically, where 
it may be shown that V = ( (N- j ) (N-k) XI + j(N-k)X2 + (N-j)kYl + jkY2) / 
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(N**2) . This can be implemented as shown in Appendix #8. In all this 
uses 6 load-from-memory, 6 multiply. 1 divide, 3 copy/add/sub 
instructions, in total 16 instructions, so again, no improvement. 

in conclusion, bilinear interpolation using 1-D LUT needs 16 
instructions to calculate v and needs 0(H) entry LUT • space . 

BILINEAR INTERPOLATION USING 2-D LUT 

As shown above, it is possible using 4 separate LUTs to lookup the 
weights to be applied to each of XI, X2, VI, Y2. using a single common 
indexing value. This index value is a function of both j and k, and each 
of these LUTs is of OW"2) size. Let these LUTs be called LUTO, LUT1, 
LUT2 , LUT3 . 

Appendix #9 shows the pseudo-code to implement this, in all this 
uses 8 load-from-memory, 5 multiply, 1 divide, 4 copy/add/sub 
instructions, in total 18 instructions. 

in conclusion, bilinear interpolation using 2-D LUT needs 18 
instructions to calculate V and needs 0(4MN«2» entry LUT-space. 

EXACT BICUBIC INTERPOLATION WITHOUT LUT 

Referring to exact 1-D cubic interpolation without LUT 24 
instructions are required to calculate value P between points wl and w2. 

Appendix #10 shows the pseudo-code to implement the technique This 
includes 24 instructions. The same code needs to be repeated for each of 
the w * Z rows similarly, to calculate points Q,R,S. So far this costs 
4*24 - 96 instructions. 

^ere is then the fifth cubic interpolation to be done in a column- 
direction, using points P.Q.R.S and tic WW as the variable. 
not need to be loaded from memory as they have ,„t S ° 4 

loads are saved on this column- interpolation. So, a further 20 
instructions are used, to give final tally of (4*24, ♦ 20 = 116 
instructions. 

in all this uses 16 load-from-memory, 30 multiply, 5 divide, 65 
copy/add/sub instructions. (If floating-point registers 
divides can be skipped at end of each row- interpolation, and the scaling 
done at end of the column interpolation, so saving a further dxvxde 
instructions, but the above allows for integer-register usage) . 



In conclusion, exact bicubic interpolation without LUT needs 116 
instructions to calculate V and needs no LUT space. 

BICUBIC INTERPOLATION USING 1-D LUT 

Referring to 1-D cubic interpolation using LUT, 16 instructions are 
required instructions to calculate value P between points wi and W2. 

Appendix #11 shows the pseudo-code to implement the technique 
including 16 instructions. The same code needs to be repeated for each of 
the w, Y, Z rows similarly, to calculate points Q,R,S. The 4 weights 
don't need to be re-loaded as they are the same, so 4 loads are saved on 
each of these 3 rows. So far this costs 16 + (3*12) = 52 instructions. 

There is then the fifth cubic- interpolation to be done in a column- 
direction, using points P,Q,R,S and t=k (Rk] as the variable. P,Q,R,S 
don't need to be loaded from memory as they have just been calculated, so 
4 loads are saved on this column- interpolation, but 4 loads are needed to 
get new weights. So, a further 12 instructions are used to give final 
tally of 16+<3*12)+12=64 instructions. 

In all this uses 24 load-f rom-memory, 20 multiply, 5 divide, 15 add 
instructions. 

(If floating-point registers are used, 4 divides can be skipped at 
end of each row- interpolation, and the scaling done at end of the column 
interpolation, so saving a further 4 divide instructions, but integer- 
register usage is allowed for) . 

In conclusion, bicubic interpolation using 1-D LUT needs 64 
instructions to calculate v and needs 0(4*N) entry LUT- space. 

BICUBIC INTERPOLATION USING 2-D LUT 

This is similar to the bilinear interpolation using a 2-D LUT 
technique, except 16 luts are needed, each with the weights to be applied 
to the 16 points wo through to Z3. 

Appendix #13 shows the technique. In all this uses 32 load-f rom- 
memory, 17 multiply, 1 divide, 16 add instructions, total 66 
instructions. Floating-point math would save the final divide, by storing 
pre -scaled weights in the LUTs. 

In conclusion, bicubic interpolation using 2-D LUT needs 66 
instructions to calculate v and needs 0(16*{N**2)) entry LUT-space. 
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HYBRID- LINEAR -BICUBIC INTERPOLATION IN GENERAL 

Referring to Figure 4, it is clear from earlier sections that 
points XI, X2, Yl, Y2 have the most importance in determining the value 
V. This is because they are the closest points to v, and also because the 
bicubic rubber- sheet passes precisely through these 4 points. To a first 
approximation, if the values at each of XI, x2, Yl, Y2 were increased by 
some small delta, then we expect the interpolated value at V to increase 
by an amount very close to delta. Thus, it can be said that alterations 
to the values at XI, X2, Yl, Y2 have a first order effect on V. 

As part of calculating V, Q and R are first calculated. Consider 
the calculation of point Q. Again. XI and X2 have a first -order effect on 
Q. points XO and X3 are used only to help calculate the slope of the 
curve at XI and at X2. Consider the calculation of slope at point XI. For 
ease of calculation and slope -continuity, the slope is typically chosen 
to be a symmetrical function, and one such function sets the slope to be 
the same as the slope of a straight line between xO and X2, i.e. slope at 

XI = (X2-X0)/2N. 

Other functions could equally well be chosen, for example an arc of 
a circle passing through XO, XI, X2 would also do, though it's more 
complicated to calculate, it can be seen that the exact choice of slope 
at XI is somewhat arbitrary, and that almost any -reasonable- choice will 
do we can never know whether the -true- original curve through xO, xl, 
X2, X3 was a cubic, a sine-wave, or some other analytic or even non- 
analytic function. 

So it can be seen that varying XO changes the slope at Xl, and has 
a second-order effect on the calculation of Q, not a first-order effect. 
The same applies to X3, YO, and Y3. 

The same argument applies to the calculation of points P and S. P 
and S each have second-order effects on v, not first-order. 

in the pure bicubic interpolation, as described above, points P and 
S are calculated using cubic interpolation, which is expensive. Consider 
point P. Points wi and W2 have a first-order effect on P, and so have a 
second-order effect on V, as point P itself only has 2nd-order effect 
on V. 

Points WO and W3 have second -order effects on P, and so have only 
3rd order effects on V. The same applies to points zO and Z3 . it is 
hardly worth wasting significant instructions calculating these 3rd-order 
effects, as 2nd-order assumptions and approximations (the slope at 
XO,xl,YO.Yl assumptions) will dwarf any such 3rd-order effects. 
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So, the present invention ignores these outermost points WO, W3, 
20, Z3, and uses linear interpolation to calculate points P and S. This 
leads to significant simplification and hence to lower implementation 
costs, both in instructions saved, and bandwidth required to load values, 
and hence to less silicon real-estate being required. 

As the significant heart of the technique is still cubic -based, 
there is no detectable loss in quality as compared to "pure" bicubic 
interpolation. 

Three implementations of the invention will be described below: 

EXACT HYBRID -LINEAR -BICUBIC INTERPOLATION WITHOUT LOT 

Referring now to Figure 6, j and k are assumed to be known, as is 
N, the inter -pixel spacing. 

j, k and N are held in registers Rj, Rk, RN respectively. 
RP and RS are two registers, just like the other R0,R1 ... 
registers. 

Wl through to 22 are not yet known, and need to be loaded from the 
image memory. 

Appendix #12 shows the pseudo-code to implement the technique. 

After getting S along the Z1-Z2 row, N*P will be in register RP, 
and N*S in register RS. So far, n*p and N*S calculations have cost 12 
instructions. 

Cubic interpolation is then applied to get Q and R, along X and Y 
rows respectively. This is similar to code described in relation to exact 
bicubic interpolation without LUT, with slightly different scaling at the 
end to calculate N*Q rather than Q. N*Q calculation has cost 24 
instructions. 

Similar code is run to calculate N*R on the Y row, costing another 
24 instructions. 

Using n*P, N*Q, N*R, N*S in registers calculated above, a column- 
direction cubic interpolation is done to calculate v. 4 load instructions 
are saved, and the final division will be division by 2(N**4). This 
interpolation of v costs 20 instructions. 

Overall cost is 12 + 24 +24 + 20 * 80 instructions, which is broken 
down into 12 load- from -memory, 22 multiply, 3 divide, 43 copy/add/sub 
instructions. 
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In conclusion, the embodiment of the invention without LUT needs 80 
instructions to calculate V and needs no LUT space. 

HYBRID -LINEAR -BICUBIC INTERPOLATION USING 1-D LUT 

This employs a combination of the code described in relation to 
both bilinear interpolation using 1-D LUT, and bicubic interpolation 
using 1-D lut. 

Appendix #14 shows the pseudo-code to implement the technique. 
Having calculated S along the zl-22 row, N*P is already in register RP, 
and N*S in register RS. n*p and N*S calculations have cost 11 
instructions. Cubic interpolation is then applied to get Q and R, along X 
and Y rows respectively. This is similar to the code described in 
relation to bicubic interpolation using 1-D LUT - again note slightly 
different scaling to calculate N*Q rather than Q. N*Q calculation has 
cost 16 instructions. 

Similar code is run to calculate N*R on the Y row. costing another 
12 instructions, as the 4 lut values are re-used, they don't need to be 
reloaded. 

using n*p, N«Q, N*R, N«s in registers calculated above, a column- 
direction cubic interpolation is used to calculate V. Again 4 load 
instructions are saved, but 4 LUT lookup- instructions need to be done. 
The final division will be division by 2(N**4) . This interpolation of v 
costs 12 instructions. 

Overall cost is 11 + 16 +12 + 12 = 51 instructions, which is broken 
down into 21 load-from-memory, 16 multiply, 3 divide, 11 copy/ add/sub 
instructions. 1 LUT is needed for the linear- interpolation, and 4 for the 
cubic stages. 

in conclusion, the embodiment of the invention using 1-D lut needs 
51 instructions to calculate v and needs 0(5*N)) entry LUT-space. 

HYBRID-LINBAR-BICUBIC INTERPOLATION USING 2-D LUT 

This is similar to bicubic interpolation using 2-D LUT. But now 
only 12 points and 12 weights need to be looked up. not 16. 

Appendix 815 shows the pseudo-code to implement the technique. This 
requires 24 load-from-memory. 13 multiply, 1 divide, 12 add instructions, 
total is 50 instructions. Floating-point math would save the final 
divide, by storing pre-scaled weights in the LUTs. 
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in conclusion, hybrid- linear-bicubic interpolation using 2-D LUT 
needs 50 instructions to calculate v and needs 0(12* (N**2)) entry LUT- 
space. 

Overall conclusion 

Figures 7 to 9 illustrate the above conclusions. 

In all cases using 1-D LUTs brings an improvement in instruction 
count, significantly so in the case of the Bicubic interpolation and in 
the case of the invention. The methods using 2-D LUTs do not have any 
significant instruction- savings over those using 1-D LUTs, and as they 
are much more expensive in space requirements (e.g. if N=256, bicubic 
using 1-D LUT takes 256 entries, against 16*256*256 « 1,048,576 entries) 
2-D LUTs are probably not worth implementing, unless their lower cost in 
multiplications or greater potential to be parallelised is significant. 

So the use of 1-D LUTs is recommended, and their storage -use of 
order N entries is very small. 

Bicubic versus Bilinear interpolation 

The fastest Bilinear uses 16 instructions, of which 6 are memory 
fetches. The fastest Bicubic uses 64 instructions, of which 24 are 
memory- fetches. Bicubic is 4x more costly than Bilinear for 
instructions, and is 4x more costly for memory- fetches . 

Invention versus Bilinear 

The fastest Bilinear uses 16 instructions, of which 6 are memory 
fetches. The best embodiment uses 51 instructions, of which 21 are 
memory- fetches. Hybrid- linear -bicubic is 3.2x more costly than Bilinear 
for instructions, and is 3 . 5x more costly for memory- fetches . 

Invention versus Bicubic interpolation 

The best embodiment uses 51 instructions, of which 21 are memory- 
fetches. The fastest Bicubic uses 64 instructions, of which 24 are 
memory- fetches. Bicubic is 1.25x more costly than the invention for 
instruction, and is 1.14x more costly for memory- fetches. 

So, the invention uses 80% of the instructions and 88% of the 
memory fetches of the traditional bicubic technique. This is a 
significant saving. The only penalty is the use of 5 rather than 4 
lookup- tables, which is insignificant. 
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It will be seen that the present invention need not be implemented 
in exactly the manner described above. For example, the linear 
interpolation could be done in the k direction, Figure 6. with the final 
cubic interpolation being done in the j direction. It will be seen, 
however, that where points are acquired in rows running in the j 
direction, eg a television signal, it is possible using the preferred 
embodiment to carry out the interpolations on the first three rows as 
they are acquired, thus reducing the time after point values are acquired 
before an interpolated result is generated. 

it will also been seen that the spacing of points on the rows need 
not be the same. This of course makes the mathematics more complex but 
does not take away from the operation of the invention which combines 
linear and cubic interpolation in an efficient yet accurate manner. 

Similarly, the points P.Q,R.S need not be located in a row running 
at 90 degrees to the w.x.Y.z rows - any transverse row will do, although 
accuracy is reduced as the P.Q.R.S row becomes more parallel to the 
w,x,Y,z rows. 

The above description has related to the implementation of the 
invention in a 2-dimensional manner, it will be seen, however, that the 
invention applies equally to systems of higher dimension. 

using either the prior art or the present invention, a 3D 
(tricubic) interpolation can be performed by reduction to a 2D problem, 
then reduction from 2D to OD (the point 'V'), or alternatively by 
reduction from 3D to ID, followed by cubic interpolation from ID (4 
points) to OD (point 'V). Similar techniques apply to yet higher 
dimensions. 

Prior art Tricubic interpolation 3D to ID then ID to OD 

using each of the four layers (each including a grid as in Figure 
3) in turn, perform 4 bicubic interpolations to reduce from 4*16 = 64 
points down to 4 new points, then perform one cubic interpolation to 
reduce from 4 new points to the final point. 64 points contribute to this 
calculation. 

Prior art Tricubic interpolation 3D to 2D then 2D to OD 

Using each of the 16 parallel rows in turn, perform 16 cubic 
interpolations to reduce from 16-4 = 64 points down to 16 new points, 
then perform one bicubic interpolation to reduce from 16 new points to 
the final point. 64 points contribute to this calculation. 
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Tricubic interpolation 3D to ID then ID to OD using the Invention 

Using the middle 4 points of each of the outer two layers, perform 
2 bilinear interpolations to reduce from 8 points to 2 new points, and 
using the two middle layers, perform 2 hybrid- linear-bicubic 
interpolations according to the invention to reduce from 2*12 = 24 points 
down to 2 new points, then perform one cubic interpolation to reduce from 
these 4 new points to the final one point. 32 points contribute to this 
calculation. 

Tricubic interpolation 3D to 2D then 2D to OD using the Invention 

Choose 4 sides of the cube which form a hollow tube, and using the 
middle 4 points of each of these 4 sides perform 8 linear interpolations, 
all going in the same direction, so that we now have 8 co-planar new 
points. Using the remaining 2 faces of the cube, choose the central 4 
points of each of these 2 faces, and also the central 4 points of the 
middle 2 layers. Using these 16 points (which form 4 rows of 4 points 
each), perform 4 cubic interpolations to calculate 4 more new points. 
These 4 new points are co-planar with the 8 points generated earlier. 
Using these 12 co-planar points, perform a 2D hybrid- linear -bicubic 
interpolation according to the invention to calculate the final point. 32 
points contribute to this calculation and so the savings of using the 
invention for either method are clear. 

it will be seen that the invention can be implemented in many forms 
and applications, such as in a digital signal processing chip, in 
graphics application software or anywhere one needs to interpolate 
information. 



18 



Appendix #1 Exact 1-D Linear Interpolation without LUT 



R2 
R2 



CODE 
Rl = Wl 
= W2 
*= Rj 
R3 = RN 
R3 -= Rj 
Rl *= R3 
Rl += R2 
Rl /= RN 



COMMENTS 

Load wl from memory 
Load W2 from memory 
R2 = W2*j 

Copy from register 
R3 = (N-j) 

Rl » (N-j)*Wl 
Rl * (W2*j + <N-j)*Wl) 



N*V 



Rl = <((W2*j + (N-j)'Wl)) / N) = V 



Appendix #2 Exact 1-D Linear Interpolation with LUT 

CODE COMMENTS 

Rl = wl Load wl from memory 

R2 = W2 Load W2 from memory 

R2 *= Rj R2 = W2*j b 
R3 = LUT [Rj 1 Lookup (N-j) in one indexed load 

Rl *= R3 Rl = (N-j)*Wl 

Rl += R2 Rl = (W2*j + <N-j)*Wl) » N*V 

■ S /-S S - U<W2*j + W-j)*HD) /N) -V 

Appendix #3 Exact 1-D Cubic Interpolation without LUT 

CODE COMMENTS 

RO = WO Load WO from memory 

Rl = wl Load Wl from memory Rl ■ d 

R2 = W2 Load W2 from memory 

R3 = W3 Load W3 from memory 

R3 -= RO R3 = (W3 - WO) 

R4 = R2 copy. R4 = W2 

R2 -= Rl R2 = (W2-W1) 
R3 -= R2 

11 - = R2 R3 = (W3-W0) - 3(W2-W1) « 2(N**3)A 

R4 RO R4 = (W2-W0) = 2*N*C 

R2 -= R3 R2 = (W2-W1) - 2(N**3)A 

we have A.B.C.D scaled by powers of N. Evaluate at t- j . 

R3 *= Rj R3 = 2(N**3)A * j 

R2 *- RN R2 - 2(N**3)B 

R3 += R2 R3 = <2(N**3)Aj + 2 (N**3)B 

11 Z V-li^^^lTiSer a instant or in sea re g ister 

R ? * = K-.iPm K)D2 N» either a constant or in some reg'r 

£ I" 3 7« 2(N**3)*Aj + 2(N**3)B)*j) * 2(N-3)C)*j) + 2<H«3)D 
R % 2 W-3) R3 i A j«3) 1 i(J«2) + C(j) ♦ D -V =Value wanted. 
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Appendix #4 Exact 1-D Cubic Interpolation with lut 

CODE COMMENTS 
RO = WO Load WO from memory 
5 Rl = Wl Load Wl from memory Rl = D 

R2 = W2 Load W2 from memory 
R3 = W3 Load W3 from memory 
R4 * LUTO [Rj] Load weight for WO from memory 

R5 = LUTl [Rj ] 
10 R6 = LUT2 [Rj ] 

R7 * LUT3 [Rj] 

RO *« R4 multiply each w by its weight 
Rl *= R5 
R2 *= R6 
15 R3 *= R7 

RO += Rl 
RO += R2 

RO += R3 RO is the sum of all 4 weighted points typically » 2(N**3)v 
RO 1= 2(N**3) Scale the result (if using integer maths). 
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Appendix #5 Exact Bilinear Interpolation without LUT 



CODE COMMENTS 

25 Rl = XI Load XI from memory 

R2 « X2 Load X2 from memory 

R2 *= Rj R2 » X2*j 

R4 = RN Copy from register 

R4 — Rj R4 » (N-j) 

30 Rl *= R4 Rl - (N-j)*Xl 

Rl R2 Rl = (X2*j + (N-j)*Xl) * N*Q 

Repeat, very similarly, to get R along Y1-Y2 row: 

35 R2 = Yl Load Yl from memory 

R3 = Y2 Load Y2 from memory 

R3 *= Rj R3 = Y2*j 

R2 *» R4 R2 » (N-j) *Y1 

R2 + = R3 R2 - (Y2*j + (N-j)'Yl) = N*R 



Interpolate along Q to R: 



R2 *» Rk R2 - NR*k 
R4 = rn Copy from register 
45 R4 -= Rk R4 = (N-k) 

Rl *- R4 Rl *= (N-k) *NQ 

Rl R2 Rl = (NR*k + (N-k) *NQ) » (N**2)V 

Rl /= (N**2) Rl = ( ( (NR*k + (N-k) *NQ) ) / (N**2)) - V 
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Appendix #6 Alternative Exact Bilinear interpolation without LUT 



CODE COMMENTS 

Rl * XI Load XI from memory 

R2 = X2 Load X2 from memory 

R3 = Yl Load Yl from memory 

R4 = y2 Load Y2 from memory 

R5 « RN Copy N 

R6 = RN Copy N 

R5 -« Rj R5 = (N-j> 

R6 - = Rk R6 = (N-k) 

Rl *= R5 Rl - (N-j)Xl 

R2 Rj R2 = jX2 

Rl += R2 Rl = (N-j)Xl + DX2 
Rl 



R6 Rl = <N-j) (N-k)Xl + j(N-k)X2 



R3 *= R5 R3 = (N-j)Yl 

R4 *= Rj R4 = jY2 

R3 R4 R3 = (N-j)Yl + JY2 

R3 *= Rk R3 « (N-j)kYl + }kY2 

R3 +• Rl 

R3 /= (N**2) R3 = V 



Appendix #7 Exact Bilinear Interpolation with LUT 

CODE COMMENTS 
Rl =* XI Load XI from memory 
K2 « X2 Load X2 from memory 

R4 ^LUTlRjl^o^kJSp in one indexed load 

Rl *= R4 Rl = (N- j) *X1 

5 +- R2 Rl = (X2*j + (N-j)*Xl) =N*Q 

Repeat, very similarly, to get R along Y1-Y2 row: 

R2 = Yl Load Yl from memory 

R3 = Y2 Load Y2 from memory 

R3 *= Rj R3 = Y2*j 

R2 R4 R2 - (N-3)*Y1 

11 + - R3 R2 = (Y2*j «• <N-j)*Yl) = N*R 

interpolate along Q to R: 

R4 = = LUT [Rk] ^Lookup*^(N-k) in one indexed load 
51 % Rl « ( ( (NR*k ♦ (N-k) *NQ) ) / (N-2) ) 



Appendix #8 Alternative Exact Bilinear Interpolation with LUT 



CODE COMMENTS 

Rl = XI Load XI from memory 

R2 = X2 Load X2 from memory 

R3 - Yl Load Yl from memory 

R4 = Y2 Load Y2 from memory 

R5 = LUT [Rj ] Lookup (N-j) in one indexed load 

R6 » LUT [RkJ Lookup (N-k) in one indexed load 

Rl *= R5 Rl = (N-j) XI 

R2 *= Rj R2 = jX2 

Rl R2 Rl m (N-j) XI + jX2 

Rl *«= R6 Rl = (N-j) (N-k) XI + j(N-k)X2 

R3 R5 R3 - (N-j)Yl 

R4 ♦» Rj R4 = jY2 

R3 += R4 R3 » (N-j)Yl + jY2 

R3 *« Rk R3 = (N-j)kYl + jkY2 

Rl += R3 Rl - (N-j) (N-k)Xl + j(N-k)X2 + (N-j)kYl + jkY2 

R3 /« (N**2) R3 - V 



Appendix #9 Bilinear Interpolation using 2-D LUT 
CODE COMMENTS 



Rj 


*- N 


Start making a 2-D index 










Rj 


+= Rk 


Rj is now a 2-D index, and has N 


**2 possible values 


Rl 


= XI 


Load XI from memory 










R2 


= X2 


Load x2 from memory 










R3 


= Yl 


Load Yl from memory 










R4 


= Y2 


Load Y2 from memory 










R5 


= LUTO [Rj 


J Lookup weight for XI 


in 


one 


indexed 


load 


R6 


« LUTl [Rj ] Lookup weight for x2 


in 


one 


indexed 


load 


R7 


= LUT2 [R j ] Lookup weight for Y2 


in 


one 


indexed 


load 


R8 


» LUT3 [Rj 


] Lookup weight for Y2 


in 


one 


indexed 


load 


Rl 


*= R5 


Rl = (N-j) (N-k) XI 










R2 


*= R6 


R2 - j ( N - j ) X2 










R3 


*= R7 


R3 = (N-j)kYl 










R4 


*« R8 


R4 « jkY2 










Rl 


+«= R2 












Rl 


+- R3 












Rl 


+= R4 


sum the 4 weighted points 










Rl 


/= (N**2) 


Rl ■ V 
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Appendix #10 Exact Bicubic Interpolation without LUT 



CODE 

R0 = wo 

Rl = Wl 
R2 = W2 
R3 = W3 
= RO 

R2 
= Rl 
= R2 



R3 
R4 
R2 
R3 



COMMENTS 

Load WO from memory 
Load wl from memory Rl 
Load W2 from memory 
Load W3 from memory 
R3 = (W3 - WO) 
copy. R4 = W2 
R2 « (W2-W1) 



R3 -= R2 



R2 
- = RO 
-= R3 
Rl 
+= RO 

we have a, 
R3 *« Rj 
R2 *= RN 



R3 
R4 
R2 
R2 
R2 



R3 
R3 
R4 



R2 
Rj 



*2) 



R3 += R4 
R3 *= Rj 
Rl *= 2(N**3 
R3 += Rl 
R3 /= 2(N**3 



R3 = 
R4 = 
R2 = 
R2 = 
R2 = 
b,C, 
R3 = 
R2 - 
R3 = 
R3 = 

R4 
R3 ° 
R3 = 
) Rl 
R3 = 
) R3 



(W3-W0) 
(W2-W0) 
(W2-W1) 
(W2-W1) 
(W2-W1) 



- 3(W2-Wl) 
« 2*N*C 

- 2(N**3)A 

- 2(N**3)A 
2(N**3)A 



= 2(N**3)A 



Wl 

. v „„ - ~ - wi 

d scaled by powers of N. 



+ WO = 2(N**2)B 
Evaluate at t=j . 



2(N**3)A 
2(N**3)B 
(2(N**3)Aj 
(2(N**3)Aj 
= 2 (N**3)C 
■ (((2(N**3)Aj 
{(((2(N**3)Aj 



+ 2(N**3)B ) 
+ 2(N**3)B ) * j 

N**2 either a constant or in some 
+ 2(N**3)B )*j) + 2(N**3)C) 
+ 2(N**3)B )*j) + 2(N**3)C) *j> 



register 



- 2(N**3)D .2(N**3) either a constant or in some reg'r 
7(< 2<N**3)Aj " 2(N"3)B)*j) ♦ 2(W*3)C)*j) + 2<N«3)D 

- Mj«3)+B(j**2)+C{j)*D -P -value wanted on 1st row. 



Repeat the above for each of the x,Y,z rows to calculate points Q,R,S. 
Then, perform a cubic interpolation using points P,Q,R,S with t-k (in Rk) 
to get final point v. 

Appendix #11 Bicubic interpolation using 1- 



D LUT 



CODE COMMENTS 

WO Load WO from memory 
wl Load Wl from memory Rl - D 
W2 Load w2 from memory 
w3 Load W3 from memory 

Load weight for WO from memory 



RO 
Rl 
R2 
R3 
R4 
R5 
R6 
R7 
RO 
Rl 
R2 
R3 



LUTO [Rj ) 
= LUT1 [Rj ] 
= LUT2 [Rj 1 

T. 1 ?? multiply each w by its weight 
*= R5 
*= R6 
*= R7 



RO += Rl 

RO +- R2 

RO += R3 RO 

RO /= 2(N**3) 



is the sum of all 4 weighted points typically 
Scale the result (if using integer maths) . 



2(N**3)V 



Reoeat the above for each of the x,Y,z rows to calculate points Q,R,S, 
« P usLg thfalready- loaded weights from the LUTs. 

Then, perform a cubic interpolation using points P,Q,R,S with t-k (in Rk) 
to get final point v. 
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Appendix #12 Exact Hybrid-Linear-BiCubic Interpolation without LUT 

CODE 
RP = Wl 
5 R2 = W2 

R2 *= Rj 
R4 = RN 
R4 -= Rj 
RP * = R4 
10 RP += R2 

n*p will not be huge, so we can assume we don't need to divide by N yet. 
Repeat, very similarly, to get S along Z1-Z2 row 

15 





RS 


= Zl 


Load 


Zl from memory 




R3 


= Z2 


Load 


Z2 from memory 




R3 


*= Rj 


R3 = 


Z2*j 




RS 


*= R4 


RS = 


(N-j) *Z1 


20 


RS 


+= R3 


RS = 


(Z2*j + (N-j)*Zl) = N*S 




Do a 


Cubic 


interpolation to get Q along x : 




RO 


= XO 


Load 


XO from memory 




Rl 


= XI 


Load 


XI from memory Rl = d 


25 


R2 


= X2 


Load 


X2 from memory 




R3 


= X3 


Load 


X3 from memory 




R3 


•= RO 


R3 = 


(X3 - XO) 




R4 


= R2 


copy. 


R4 = X2 




R2 


-« Rl 


R2 = 


(X2-X1) 


30 


R3 


- = R2 








R3 


- = R2 








R3 


-» R2 


R3 = 


(X3-X0) - 3(X2-Xl) = 2 




R4 


-= RO 


R4 - 


(X2-X0) ~ 2*N*C 




R2 


-= R3 


R2 = 


(X2-X1) - 2(N**3)A 


35 


R2 


•* Rl 


R2 = 


(X2-X1) - 2(N**3)A - XI 




R2 


+= RO 


R2 = 


(X2-X1) - 2{N**3)A - XI 



We have A,B,c,D scaled by powers of N. Evaluate at t=j . 
R3 *= Rj R3 - 2(N**3)A * j 
R2 *= RN R2 = 2(N**3)B 
40 R3 += R2 R3 - (2(N**3)Aj + 2(N**3)B ) ~*~ 

R3 *= Rj R3 * (2(N**3)Aj + 2{N**3)B ) * j 

R4 *= (N**2) R4 = 2(N**3)C .N**2 either a constant or in some register 
R3 +■ R4 R3 = (((2(N**3)Aj + 2(N**3)B ) * j ) + 2<N**3)C) 
R3 *■ Rj R3 = ((((2(N**3)Aj + 2(N**3)B )*j) + 2 (M**3) C)' "♦Jj"™" 
45 Rl *- 2(N**3) Rl = 2(N**3)D .2<N**3) either a constant or in some reg'r 

R3 += Rl R3 = (<((2<N**3)Aj + 2(N**3)B)*j) + 2(N**3)C)*jj + 2<N**3)D 
R3 /= 2(N**2) R3 = (A<j**3) ♦ B(j**2) + C(j) + D) *N =N*Q . 

Repeat the above cubic code for the Y row to calculate point R. 

50 - . ._ 

Then, perform a cubic interpolation using points P,Q,R,S with t=k (in Rk) 
to get final point V. 



COMMENTS 

Load wl from memory 
Load W2 from memory 
R2 = W2*j 

Copy from register 
R4 = (N-j) 
RP = (N-j)*Wl 

RP = (W2*j + (N-j)*Wl) » N*P 
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Appendix #13 Bicubic Interpolation Using 2-D LUT 

CODE COMMENTS 
Rj *«= N Start making a 2-D index ... 
5 Rj += Rk Rj is now a 2-D index, and has N**2 possible values. 

RO = WO Load WO from memory 

Rl = wl Load Wl from memory 

R2 = W2 Load w2 from memory 

10 R3 = W3 Load W3 from memory 

R4 = XO Load XO 

etc up to Load X,Y,Z values ... 

R14 «= Z2 n 

R15 « Z3 16 Loads to get WO - Z3 

15 R i6 = LUTO [Rj ] Lookup weight for wo in one indexed load 

R17 - LUTllRj] Lookup weight for Wl in one indexed load 
... etc up to . . .more loads to get weights 

R31 = LUTl5[Rj] Lookup weight for z3 in one indexed load 

20 . _ fc 

RO *= R16 RO = WO * its weight 
Rl *= R17 Rl « Wl * its weight 

etc up to . . 16 rapys to multiply all weights by values 
R15 R31 R15 = Z3 * its weight 

RO += Rl sum the 16 weighted points 

pQ += R2 

etc up to . . 15 adds to sum 16 weighted values 

RO += R14 1 , 

RO +- R15 all 16 weighted points summed. 

RO /= 2(N**3) RO = V 
Appendix #14 Hybrid- Linear -BiCubic interpolation using 1-D LUT 

CODE COMMENTS 

rp B wl Load Wl from memory 
R2 = W2 Load W2 from memory 
r2 *= Rj R2 = W2*j 

R4 = LUT [Rj 1 Lookup (N-j) in one indexed load 
RP R4 RP = (N-j)*wi 
R p += R2 RP - <W2*j + (N-j)*Wl) = N*P 
Repeal very similarly, to get S along Z1-Z2 row, so we have 
RS = Zl Load Zl from memory 
R3 = z2 Load Z2 from memory 
R3 *» Rj R3 = Z2*j 
RS *= R4 RS = (N-j)*Zl 
RS += R3 RS = (Z2*j + (N-j) •ZD - N*S 
Do a Cubic interpolation to get Q along X row 
RO = XO Load XO from memory 
Rl = XI Load XI from memory Rl = d 
R2 « X2 Load X2 from memory 
R3 = X3 Load X3 from memory 

R4 m LUTO [Rj ) Load weight for XO from memory 
55 R5 = LUTltRjl 

R6 * LUT 2 [Rj] 

R7 = LUT 3 [Rj] ^ . . . 

RO *« R4 multiply each X by its weight 

Rl *= R5 
60 R2 *» R6 

R3 *= R7 
RO += Rl 

5! Z R3 RO is the sum of all 4 weighted points typically = 2(N**3)v 
65 Repea^thf ioove «Sc'c22 for the Y row to calculate point R. 

Then, perform a cubic interpolation using points P,Q,R,S with t=k (in Rk) 
to get final point v. 



25 



30 



35 



40 



45 



50 



25 



Appendix #15 Hybrid-Linear-BiCubic Interpolation using 2-D LUT 



CODE 
Rj *= N 
Rj += Rk 



RO .= 
Rl = 
R2 = 
R3 = 
R4 = 
R5 = 
R6 = 
R7 = 
R8 = 
R9 = 
RIO 
Rll 



Wl 
W2 
XO 
XI 
X2 
X3 
YO 
Yl 
Y2 
' Y3 
= Zl 

« Z2 



COMMENTS 

Start making a 2-D index ... 

Rj is now a 2-D index, and has N**2 possible values. 



Load Wl 
Load w2 
Load XO 
Load XI 
Load X2 
Load X3 
Load YO 
Load Yl 
Load Y2 
Load Y3 
Load Zl 
Load Z2 



from memory 



R12 = lutO [Rj ] Lookup weight for Wl in one indexed load 

R13 = LUT1 [Rj ] Lookup weight for W2 in one indexed load 
. . . etc up to . . . 

R23 « LUT15[Rj] Lookup weight for Z2 in one indexed load 



RO *= R12 RO ■ Wl * its weight 
Rl *= R13 Rl = W2 * its weight 
R2 *= R14 R2 = XO * its weight 

.. etc up to 12 mpys to multiply all weights by values 
Rll *= R23 Rll =* Z2 * its weight 

RO += Rl sum the 12 weighted points 
RO +* R2 

. . etc up to . . 12 adds to sum 12 weighted values 
RO += RIO 

RO += Rll all 12 weighted points summed. 
RO /= 2(N**3) RO = V 
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CLAIMS 



1. A method of determining an interpolated value for a point lying 
between second and third rows of an array comprising at least four rows 
of spatially related points each having respective values, said method 
comprising the steps of: 

linearly interpolating between a first two points on said first row 
to determine a first value at a first location; 

cubically interpolating between a second two points on said second 
row in accordance with the rate of change of values at said second two 
points on said second row to determine a second value at a second 
location; 

cubically interpolating between a third two points on said third 
row in accordance with the rate of change of values at said third two 
points on said third row to determine a third value at a third location; 

linearly interpolating between a fourth two points on said fourth 
row to determine a fourth value at a fourth location; said first, second, 
third and fourth locations forming a transverse row intersecting said 
point; and 

cubically interpolating between said second and third locations in 
accordance with the rate of change of values at said second and third 
locations determined according to said first and fourth values to 
determine said interpolated value. 

2 a method according to claim 1 wherein said linear interpolating 
steps employ a look-up table which provides weightings for use in said 
interpolating steps. 

3 A method according to claim 1 or 2 wherein said cubic interpolating 
steps employ a look-up table which provides weightings for use in said 
second and third row interpolating steps. 

4 An interpolation component for determining an interpolated value 
for a point lying between second and third rows of an array comprising at 
least four rows of spatially related points each having respective 
values, said component comprising: 

means instantiable to linearly interpolate between a first two 
points on said first row to determine a first value at a first location; 
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means instantiable to cubically interpolate between a second two 
points on said second row in accordance with the rate of change of values 
1 at said second two points on said second row to determine a second value 
at a second location; 

means instantiable to cubically interpolate between a third two 
points on said third row in accordance with the rate of change of values 
at said third two points on said third row to determine a third value at 
a third location; 

means instantiable to linearly interpolate between a fourth two 
points on said fourth row to determine a fourth value at a fourth 
location; said first, second, third x and fourth locations forming a 
transverse row intersecting said point; and 

means instantiable to cubically interpolate between said second and 
third locations in accordance with the rate of change of values at said 
second and third locations determined according to said first and fourth 
values to determine said interpolated value. 

5- A digital signal processor including the interpolation component of 
claim 4. 

6. A graphics application including the interpolation component of 
claim 4 . 

7. A computer program product comprising computer program code stored 
on a computer readable storage medium for, when executed on a computer, 
determining an interpolated value for a point, the program code 
comprising an interpolation component as claimed in claim 4. 
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